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It is known that point searching in basic semialgebraic sets and the search for 

globally minimal points in polynomial optimization tasks can be carried out using 

{s d)^"'' arithmetic operations, where n and s are the numbers of variables and 

^ N . constraints and d is the maximal degree of the polynomials involved. 

ry\ ', We associate to each of these problems an intrinsic system degree which becomes 

jy^ I in worst case of order {nd)^""' and which measures the intrinsic complexity of 

O ■ the task under consideration. 

We design non-uniformly deterministic or uniformly probabilistic algorithms of 
. intrinsic, quasi-polynomial complexity which solve these problems. 
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1 Introduction 

We develop uniform bounded error probabilistic and non-uniform deterministic algo- 
rithms of intrinsic, quasi-polynomial complexity for the point searching problem in 
basic semialgebraic sets and for the search of isolated local and global minimal points 
in polynomial optimization. The semialgebraic sets and optimization problems have to 
satisfy certain well motivated geometric restrictions which allow to associate with them 
an intrinsic system degree (see Section I3.2p which controls the complexity of our algo- 
rithms and constitute the core of their intrinsic character. The algorithms we are going 
to design will become then polynomial in the length of the extrinsic description of the 
problem under consideration and its system degree (we take only arithmetic operations 
and comparisons in Q into account at unit costs). The idea is that the system degree 
constitutes a geometric invariant which measures the intrinsic "complexity" of the con- 
crete problem under consideration (not of all problems like a worst case complexity). 
In worst case the sequential time complexity will be of order (^^ (n d)'^^"'^ or (nd)^^"'^ , 
where n is the number of variables and d the maximal degree of the polynomials 
occurring in the problem description, s their number and 1 < p < n the maximal 
codimension of the real varieties given by the active constraints. We shall suppose that 
these polynomials are represented as outputs of an essentially division-free arithmetic 
circuit in Q[Xi, . . . ,Xn] of size L (here, we mean by essentially division-free that 
only divisions by rational numbers are allowed). The (sequential) complexity of our 
algorithms is then of order Li^^^n^^^^d'^^^^S'^ (or L(nd)'^^^^S'^), where 6 is the intrinsic 
system degree which in worst case becomes of order {srid)^^"^^ (or (nd)'^^'^^). We call 
this type of complexity bounds intrinsic and quasi-polynomial. 

For the problem of deciding the consistency of a given set of inequality constraints 
and to find, in case the answer is positive, a real algebraic sample point for each con- 
nected component of the corresponding semialgebraic set sequential time bounds of 
order (sd)^^''^ , (srf)^("') and (s(i)'^("') are exhibited in [211 [El EH EZl ESI [71 129] and 
in the book [8]. Such bounds can also be derived from efficient quantifier elimination 
algorithms over the reals ([271 EH [71 [8]). Since two alternating blocks of quantifiers be- 
come involved, one would expect at first glance that only a (sc?)*^'-" ^ time complexity 
bound could be deduced from efficient real quantifier elimination for polynomial opti- 
mization problems. But, at least for global optimization, one can do much better with 
an s2n+i^o(n) sequential time bound (see [8], Algorithm 14.46, the result is essentially 
due to [35]). For particular global polynomial optimization problems the constant hid- 
den in this bound can be made precise and the algorithms become implementable (see 
[221 [Ml [ST] and [21] )• Accurate estimations for the minima are contained in [2U] . 

Nevertheless, this article does not focus on the improvement of known worst case 
complexity bounds in optimization theory. Our aim is to exhibit large classes of point 
searching problems in semialgebraic sets and polynomial optimization problems where 
it makes sense to speak about intrinsic complexity of solution algorithms. This is the 
reason why we put the accent on geometrical aspects of these problems. The algorithms 
become then borrowed from [20] (see also [ISl [IH1[2S1 [13 [H]) and, in particular, from 
[21 [3] (or alternatively from |36j). 



It is not straightforward but possible to transfer our algorithms and complexity results 
to the bit model. The main task of such a transfer concerns the algorithms from [21 [3] 
we use and this exceeds the scope of this paper. 

1.1 Notions and Notations 

We shall freely use standard notions, results and notations from algebraic and semial- 
gebraic geometry, commutative algebra and algebraic complexity theory which can be 
found e.g. in the books [321 EH E31 [ID] • 

Let Q , M and C be the fields of the rational, real and complex numbers, respectively, 
let Xi, . . . , Xn be indeterminates over C and let Fi, . . . , Fp be a regular sequence 
of polynomials in ]R[Xi, . . . ,Xn] defining a closed, R -definable subvariety S of the 
n -dimensional complex afiine space C"^ . Thus S" is a nonempty equidimensional afiine 
variety of dimension n—p, i.e., each irreducible component of S is of dimension n — p. 
We denote by S'k := S'nM" the real trace of the complex variety S . We shall use also 
the following notations: 

{Fi = 0, . . . , Fp = 0} := 5 and {Fi = 0, . . . , Fp = 0}m := S^. 

For a given polynomial G G M[Xi, . . . ,X„] we denote by Sg and (5*15)0 the sets of 
points of S and S'm at which G does not vanish. 

We call the regular sequence Fi, . . . , Fp reduced if for any index 1 < k < p the ideal 
(Fi, . . . , Ffc) is radical. A point x of C" is called (Fi, . . . , Fp) -regular if the Jacobian 

J(Fi, . . . , Fp) := -7^^- has maximal rank p at x . Observe, that for each reduced 



i<j<p 

l<k<n 



regular sequence Fi, . . . , Fp defining the variety S , the locus of (Fi, . . . , Fp) -regular 
points of S is the same. In this case we call an (Fi, . . . , Fp) -regular point of S simply 
regular (or smooth) or we say that S is regular (or smooth) at x . The variety S is 
called (Fi, . . . , Fp) -regular or smooth if S is (Fi, . . . , Fp) -regular at any of its points. 

For a given afiine variety V denote by C[\^] its coordinate ring. We call the elements 
of C[V] coordinate functions. For g G C[V] we denote by C[y]g the localization of 
C[V] at g. Observe that C[y]g is isomorphic to C[V^] . 

Suppose for the moment that 1/ is a closed subvariety of C" . For V irreducible we 
define its degree deg V as the maximal number of points we can obtain by cutting V 
with finitely many affine hyperplanes of C" such that the intersection is finite. Observe 
that this maximum is reached when we intersect V with dimension of V many generic 
afiine hyperplanes of C" . In case that V is not irreducible let V = CiU ■ ■ ■ U Cg 
the decomposition of V into irreducible components. We define the degree of V as 
degy :=Ei<j<sdegCj. 

With this definition we can state the so-called Bezout Inequality: 
Let V and W be closed subvarieties of C" . Then we have 

deg(y r\W)<degV ■ deg W. 



If y is a hypersurface of C" then its degree equals the degree of its minimal equation. 
The degree of a point of C" is just one. For more details we refer to [25l [T6l HO]. 

Let 1 < i < n — p and let a := [a^ i] i<k<n-p-i+i be a real ((n — p — i + 1) x (n + 1) - 

' 0<i<n 

i+1 o) 7^ and suppose that \ak i] i<k<n-p-i+i 



) 0'n—p—i-\ 



has maximal 



matrix with (aio, 

rank n — p — i + 1 . 

The dual ith polar variety of S associated with the matrix a is defined as closure of 

the locus of the (Fi, . . . , Fp) -regular points of S where all (n — z + 1) -minors of the 

polynomial ((n — i + 1) x n) -matrix 



dFi 

dFp 
dXi 

11,1 — "1,0^1 



dFi 
dFj, 

axn 

0-1,71 — 0,lfiXr, 



,0-^1 



-f — i + l,n 



-p — i + 1 



,0^n. 



vanish. 



Observe that this definition of dual polar varieties may be extended to the case that 
there is given a Zariski open subset O of C" such that the equations Fi = 0, . . . , Fp = 
intersect transversally at any of their common solutions in O and that S is now the 
locally closed subvariety of C" given by 

S- := {Fi = 0, . . . , Fp = 0} n O. 

In [2] and [3] we have introduced the notion of dual polar variety of S and motivated 
by geometric arguments the calculatory definition above of these objects. Moreover, we 
have shown that, for a real {{n — p — i + 1) x (?i + 1)) -matrix a = [ak,i] i<k<n-p-i+i with 

0<l<n 

[ofc i] i<fc<„-p-i+i generic, the zth polar varieties is either empty or of pure codimension 

' l<l<n 

i in S . Further, we have shown that this polar variety is normal and Cohen-Macaulay 
(but not necessarily smooth) at any of their (Fi, . . . , Fp) -regular points (see |1], Corol- 
lary 2 and Section 3.1). This motivates the consideration of the so-called generic dual 
polar varieties associated with real {{n — p — i + l)x[n + 1)) -matrices a which are 
generic in the above sense, as invariants of the variety 5* (independently of the given 
equation system Fi = 0, . . . , Fp = ). 

For our use of the word "generic" we refer to [1], Definition 1. 

In case that S^ is smooth and a is a real {{n — p — i + 1) x (n+ 1)) -matrix, as 
before, the i th dual polar variety associated with a contains at least one point of each 
connected component of ^jr and is therefore not empty (see [2] and [3], Proposition 2). 



2 Inequalities 

Let Fi, . . . , Fg G M[Xi, . . . , X„] and I < p < min{s, n} . We assume that for any 
1 < j'l < ■ ■ ■ < jfc < s, 1 < k < p the following condition is satisfied. 



Condition A 

The polynomials Fj-^, . . . ,Fj^ generate in M.[Xi, . . . ,Xn] the trivial ideal or form a 
reduced regular sequence and the semialgebraic set {Fj-^ = 0, . . . , F^-^ = 0}r is empty or 
{Fj^, . . . , FjJ -regular. Moreover, any p+1 polynomials of Fi, . . . ,Fs have no common 
real zero. 

For El, . . . ,es E { — 1, 1} let 

{signFi = El,. . .,signFs = ej := {x G M" | signFi(a;) = ei, . . . ,signFs(x) = Eg}. 
From now on we shall suppose without loss of generality Ei = ■ ■ ■ = Eg = 1 and write 

{Fi > 0,. . .,Fs> 0} instead of {signFi(a;) = 1, . . . ,sign Fs{x) = 1}. 
Let us fix a maximal index set I < ji < ■ ■ ■ < jk ^ s, 1 < k < p with 



{F,, = 0, . . . , F,, = 0}m n {Fi > 0, . . . , F, > 0} 7^ 



Lemma 1 



{Fj, = 0, . . . , F,, = 0}m n {Fi > 0, . . . , F, > 0} = 

= {Fj, = 0, . . . , F,, = 0}m n {F,- > 0, 1 < J < s, jV Ji, • • • , jV 3k}. 

Proof. We show first the inclusion of the left hand side of the set equation in the 
right hand side. For this purpose, let x be an arbitrary point of {Fj-^ = 0, . . . , Fj^ = 
0}iK n {Fi > 0, . . . , Fg > 0} . Suppose that there exists an index I < j < s, j j^ 
ji,...,j 7^ jk with Fj{x) = 0. Then x belongs to {Fj^ = 0, ...,Fj^ = 0,Fj = 
0}r n {Fi > 0, . . . ,Fs > 0} which by the maximal choice oi 1 < ji < ■ ■ ■ < jk < s is 
empty. Therefore, we have Fj{x) ^ for any index 1 < j < s, j y^ ji, . . . ,j ^^ jk- 
Since x belongs to {Fi > 0, . . . , F^ > 0} , we have Fj{x) > . 

We are now going to show the inverse inclusion. Consider an arbitrary point x G 
{Fj, = 0, . . . , F,-, = 0}m n {Fj >0, 1<J <s, j ^n,...,j ^ Jk} and let U be an 
arbitrary neighborhood of x in M" . Without loss of generality we may assume that 
U is contained in {Fj > 0, 1 < j < s, j j^ ji, . . . ,j j^ jk} ■ Since by Condition A 
the point x is contained in the (F^-^, . . . , Fj^. ) -regular set {Fj^ = 0, . . . , F,,, = 0}k , the 
polynomial map from M" into M" given by (F,-^, . . . ,Fj^) is a submersion at x and 
therefore there exists a point y E U with Fj^iij) > 0, . . . , Fj^iij) > . Because U was 
an arbitrary neighborhood of x, we conclude that x belongs to {Fj^ = 0, . . . , F,-^ = 
0}m n {Fi > 0, . . . , F, > 0} . a 

Corollary 2 

Let C be a connected component of {Fj-^ = 0, . . . , Fj^, = 0}^ with Cfl {Fi > 0, . . . , F^ 
>O}^0. Then 

C C {F,- > 0, 1 < i < s, i ^ ii, . . . ,i ^ 3k}. 



Proof. Lemma [T] implies that the set 



{F,, = 0, . . . , F,, = 0}m n {Fi > 0, . . . , F, > 0} 

is open and closed in {Fj^ = 0, . . . , Fj^ = 0}k . Therefore this set is the union of all the 
connected components of {F,-^ = 0, . . . , Fj^. = 0}r which have a nonempty intersection 
with {Fi > 0, . . . , F, > 0} . This imphes Corollary El D 

2.1 Converting nonstrict inequalities into strict ones 

Let I < k < p, 1 < ji < ■ • ■ < ife < s and let x = {xi, . . . ,Xn) be a point of 
{Fj-^ = 0, . . . , Fji_ = 0}]R n {Fi > 0, . . . ,Fs > 0} . Thus x satisfies the system of 
nonstrict inequalities 

F,-,(x) > 0, . . . , F,„(x) > 0, F,(x) > 0, 1 < J < s, jV Ji, • • • , jV Jk 

Starting from x we wish to construct a point y G M" which satisfies the strict inequal- 
ities 

F,iy)>0, l<j<s. 

From Condition A we conclude that the Jacobian J{Fj^, . . . , F,^.) has full rank k ai x . 
Therefore we may efficiently find a vector /j, = (//i, . . . , /i„) G M."^ such that the entries 
of J{Fj^, . . . , Fj^){x)ii^ are all positive (here, jjF denotes the transposed vector of /x). 

Let y be a new indeterminate and for 1 < j < s let Gj := Fj{fiiY + xi, . . . , /U„y + 
Xn) ■ Observe, that the univariate polynomial Gj satisfies the equation -^(0) = 
J2i<i<n dx'i^)f^i ■ In particular, the entries of 

(^(0),...,^(0)) = J(F,„...,F,J(x)/.^ 

are all positive. Let c > be the smallest positive zero of HkiXs ^j (if there exists 
none, c may be any positive real number). Then one verifies immediately that z := 
a; + |/i satisfies for any index 1 < j < s the condition Fj{z) = Gj(|) > 0. 

2.2 Finding sample points for all consistent sign conditions 

Let {si, . . . ,es) G { — 1, 0, 1}* . The polynomial inequality system sign Fi = ei, . . . , sign 
Fg = Eg is called a sign condition on Fi, . . . , F, which we say to be consistent if there 
exists a point x G M" satisfying it. In case (ei, . . . , e^) G { — 1, l}** we call the sign 
condition strict, otherwise nonstrict . A real algebraic point of M" which is supposed 
to be encoded "a la Thom" [13] and to satisfy the sign condition is called a sample 
point of it . 

Let Fi, . . . ,Fs be given as outputs of an essentially division-free arithmetic circuit /3 
in Q[Xi, . . . , Xn] having size L and non-scalar depth i (hence, Fi, . . . , F^ belong to 
Ti, . . . , Xn] ). Let d > 2 be an upper bound of deg Fi, . . . , deg Fg . For 1 < k < p 

6 



and 1 < ji < ■ ■ ■ < ifc < s let 5ji,...,j,, be the maximal degree of {Fj^ = 0, . . . , Fj^, = 0} 
and all generic dual polar varieties of this variety. Let finally 

6 := max{5j^^...j^ \ I < ji < ■ ■ ■ < jk < s, 1 <k <p}. 

We call 6 the degree of the sample point finding problem for all consistent sign con- 
ditions of Fi, . . . ,Fs . From the Bezout Inequality we deduce 6 < p^~Pcir- . Using 
the algorithms from [21 |3] (see also [5] for technical details) we construct for each 
1 < k < p and 1 < ji < ■ ■ ■ < jfc < s real algebraic sample points for each con- 
nected component of {Fj^ = 0, . . . , F^,. = 0}m. Then we evaluate the signs of all 
Fj, 1 < J < s j 7^ ji, . . . , j 7^ j'fc on these sample points. In this way we obtain sample 
points for all nonstrict consistent sign conditions on Fi, . . . ,Fs . As far as only sample 
points for the strict sign conditions on Fi, . . . , Fg are required, we limit our attention 
to sample points of the connected component of {Fj^ = 0, . . . , Fj^ = 0}^ where the 
signs of aAl Fj,l < j < s j y^ ji, ■ ■ ■ ,j j^ jk are all strict. Let x be such a sample 
point with sign Fj{x) = Sj , and £j G { — 1, 1} for 1 < j < s with j 7^ ji, ■ ■ ■ ,j 7^ jk- 
Then, following Subsection 12.11 for any Ej-^, . . . ,ej^. G { — 1, 1} we may convert x into 
a real algebraic sample point of the strict sign conditions signFi = £1, . . . ,signFs = 
Eg- The whole procedure can be realized in time L(^) nP^^'' dP^^^ 5"^ (here arithmetic 
operations and comparisons in Q are taken into account at unit costs). We have 
therefore shown the following statement which constitutes a simplified variant of [21], 
Theorem 5. 

Theorem 3 

Let n, d, p, s, L, i, 6 G N with 1 < p < n be arbitrary and let Fi, . . . , Fg G 
Q[Xi, . . . , Xn] be polynomials of degree at most d satisfying Condition A and having 
sample point finding degree at most 6 . Suppose that Fi, . . . , Fg are given as outputs 
of an essentially division-free circuit /3 in Q[Xi, . . . ,Xn] of size L. 

There exists a uniform bounded error probabilistic algorithm A over Q which com- 
putes from the input (3 in time L(*) rfi^v) d^^^^ S"^ = s'^^^^nd)'^^"'^ real algebraic sample 
points for each consistent sign condition on Fi, . . . ,Fs (here, algebraic operations and 
comparisons in Q are taken into account at unit costs). 

For any n, d, p, s, L, £, 6 ^N with 1 < p < n the probabilistic algorithm A may be 



realized by an algebraic computation tree of depth L ('^) n'~'(p) d'^(i) S'^ = s'~'^\ndY 



that depends on parameters which may be chosen randomly. 



Observe that in the elimination process which leads to Theorem [3l the intermediate 
polynomials may become of degree nd6 . Thus 6 is not an upper bound for them. 

3 Optimization 

We associate with a polynomial optimization problem with smooth equality constraints 
certain natural geometric conditions and an intrinsic invariant that controls the com- 
plexity of the algorithm which we are going to develop in order to solve this problem. 
Our approach has some features in common with that of 



3.1 Geometric considerations 

Let be given polynomials G,Fi, . . . , Fp E M[Xi, . . . , Xn], I < p < n, and suppose that 
Fi, . . . , Fp form a reduced regular sequence and that the real trace V^ of 

y := {x e C" I Fi(x) = 0, . . . , Fp{x) = 0} 



is nonempty and smooth. 

For I < ji < ■ ■ ■ < jp < n denote by Aj^_...jp the p-minor of the Jacobian J(Fi, . . . , Fp 

given by the columns numbered ji, . . . ,jp and by 

]\Ajl,--;jp) ]\Ajl,--;jp) 

IVl^ , . . . , IVl^-p 

the {p + 1) -minors of ((p + 1) x n) -matrix 

■j(Fi,...,F,) 

dG dG 

. dXi" ' dX„ 

given by the columns numbered ji, . . . , jp to which we add, one by one, the columns 

numbered by the indices belonging to the set {1, . . . , r?,} \ {ji, . . . ,jp} . 

Let a; be a local minimal point of G on V^ • Then the Karush-Kuhn- Tucker conditions 

dG 



imply that rk 



dXi 



■.Fp) 

dG 
dX„ 



(x) < p holds. Let 



W ■.= {xeV\rk 



J{F,,...,Fp) 



dG 
L dXi 



dG 
dX„ J 



(x) < p} 



and let 1 < ji < ■ ■ ■ < jp < n. Then we deduce from the Exchange Lemma in [T] that 



Wa 



■31,---,3p 



Va 



Ji,---,]p 



n{M\ 



(iivjp). 



X) 



0, 



M 



in, -dp) 



{x) = 0} 



holds. 

From now on we assume that G and Fi, . . . ,Fp satisfy the following three conditions 

for any 1 < ji < ■ ■ ■ < jp < n and any 1 < k < n — p. 

Let Dk denote the union of all irreducible components of W of dimension strictly 
larger than n — p — k . Observe that Dn-p+i is also well defined. 

Condition B 

Any point of 



(Mr 



^J\,---.3p 

r(jivjp) 



n {m[''-''\x) = 0, . . . , m]^''-''\x) 



0}r \(Dfc)i 



m'^^^'-'^"^) -regular. 



IS {F^, . . . , Fp, Ml' 

Condition B implies that the real trace of V^j^,...,, njAfl^''''^"^ 
is smooth. Let Ci, . . . , C^ be the irreducible components of 



M, 



Ul,--;jp) 



0} 



(y^n....P n {^1 



ijl,--;jp) 



0,...,M^- 
8 



(iivjp) 



0})\D, 



Condition C 

For any 1 < j < s the semialgebraic set 

(C,-)r \ (Ci U ■ ■ ■ U Q_i U C,+i U ■ ■ ■ U Cs)r 

is dense in (Cj)m and {Cj)r is nonempty. 

Condition C implies tliat each irreducible component of W contains a real point, 
that (Ci)]R, . . . , (Cs)r are the irreducible components of the real trace of (Va^^ . ,,■ n 
^M(n,-,jp) = 0, . . . , MJ^^^'-'^"^ = 0}) \ Dk and that Ci, . . . , C, can be obtained from 
the complex interpretation of the defining equations of (Ci)]r, . . . , {Cs)m. ■ 

Condition D 

The polynomial Ajj^...jp does not vanish identically on any irreducible component of 

V n {mJ^'''-'^"^ = 0, . . . , mI^''-'^"'^ = 0} . 

We need Condition D only for the proof of Lemma [6] below. We may therefore replace 
Condition D by the statement of Lemma O Condition D is the less intuitive one of our 
conditions. We need it to make the algorithm in Section 13.21 work correctly. 

Lemma 4 

Let notations be as before. For any 1 < j < s we have dim Cj = dim(Cj)]K = n—p — k . 

Proof. Since (Cj)ir is nonempty there exists by Condition C an open semialgebraic 

subset t/ of M'" , disjoint from (DfcUCiU- ■ -UQ^iUCj+iU- • ■UCs)r , with f/n(Cj)iR 7^ . 
This implies 

u n {Vr)a,^_,^ n {m['"-''''\x) = 0, . . . , Mi''-''\x) = o}^ = un (C,)m. 

From U fl Dk = and Condition B we infer now that any point of the semialge- 



(iivjp) ]\/r{Ji,---Jp)\ 



braic set U n {Cj)^ is (Fi, . . . , Fp, Ml^''-'^"' , . . . , Ml^''-'^"') -regular. Hence, dim C. 



dim(Cj)iR = n — p — /c. □ 

Lemma 5 

Let C be an irreducible component of W . Then G takes a constant real value on C . 

Proof. Since by assumption V^ = {Fi = 0, . . . , Fp = OJk is smooth and since 
Ck is nonempty by Condition C, there exist indices I < ji < ■ ■ ■ < jp < n with 
(C*ir)a 7^ 0- Furthermore, there exists an index 1 < k < n — p with dimC = 

n — p — k . 

From now on let notations be as before. Without loss of generality we may assume that 
Caj^,...j \Dk = Ci holds. By Condition C the semialgebraic set (Ci)ir\(C2U- ■ ■UCs)r 
is nonempty. Let x be an arbitrary point of this set. As we have seen in the proof of 
Lemma HI there exists an open semialgebraic neighborhood f/ of z in W^ with 

u n (M8)a,,....„, n {mJ^'^'-*) = 0, . . . , Mi'"-''"^ = o}k = f/ n (Ci)m 

and Un{Dk)R = 0- Condition B implies now that Un{Ci)R is a smooth semialgebraic 
manifold which we may suppose to be connected by continuously differentiable paths. 

9 



Let y be an arbitrary point of ^ fl (Ci)m and let r be a continuously differentiable 
path in t/ n (Ci)]r that connects x with y . We may suppose that r can be extended 
to a suitable open neighborhood of [0, 1] in M and that r(0) = x and r(l) = y 
holds. Observe that r([0, 1]) is contained in V^. The path r depends on a param- 
eter T defined in the given neighborhood of [0,1]. Let t{T) := (ri(T), . . . , r„,(T)) . 
Since r([0,l]) is contained in Vm. the vector (^(t), . . . , ^(t)) belongs to the ker- 
nel of J((Fi,...,Fp))(r(t)) for any t G [0,1]. On the other hand, Ci C C C W 
implies that (^^(r(t)), . . . , ^-(r(f)) is linearly dependent on the full rank matrix 

J(Fi,...,Fp)(r(/)). Therefore we have ^(t) = ELi ffi(^W)§W = for any 
t G [0, 1] . Hence, Got is constant on [0, 1] . Consequently, we have G{x) = G{t{0)) = 
(j(r(l)) = G{y) . From the arbitrary choice of x and y in U H {Ci)m. we infer that G 
takes on U (1 (Ci)ir a constant value. Thus the restriction of G to the semialgebraic 
set (Ci)iR \ (C2 U ■ ■ ■ U Cs)r is locally constant and takes therefore only finitely many 
values in M . Condition C implies now that the same is true for the restriction of G to 
(C'i)r- 

By Conditions B and C there exists an (Fi, . . . , Fp, m[^''-'^''\ . . . , MJ^^^'-'^"^) -regular 
point X = {xi, . . . ,Xn) of (Ci)k \ (C2 U ■■■ U Cs)r. Hence, there exists an open 
semialgebraic neighborhood ^ of x in M" and n — p — k parameters ^1, . . . ,^n-p-k of 
UnGi such that the restriction of G to Un {Ci)r can be developed into a convergent 
power series P{Ci ~ ^^i, . . . , in-p-k — Xn-p-k) around (xi, . . . , Xn-p-k) ■ Since G is 
locally constant on f/ fl (Ci)m we conclude that P{C,i — Xi,..., ^n-p-k — ^n-p-k) equals 
its constant term, say 6 G M . 

On the other hand, there exists an open neighborhood O of x in C" such that the 
restriction of G to O fl Ci can be developed into a convergent power series in ^1 — 
Xi, . . . , Cn-p-k — Xn~p~k ■ This powcr series must necessarily be P{Ci ~ 2;i, . . . , in-p-k — 
Xn-p~k) ■ Thus G takes on O fl Ci only the real value b. Suppose that G takes on 
Ci a value different from h. Then {Cijc^b is nonempty and therefore (by [32], Ch. I, 
Section 10, Corollary 1) dense in the Euclidean topology of Ci . In particular, there 
exists a point y G O fl Ci with G{y) 7^ h. This contradiction implies Lemma [51 □ 

Let 1 < k < n — p. By Lemma [5] the polynomial G takes on Dk only finitely 
many values which are all real. Denote by B^ the set of these values and let Rk '■= 
YlbesX^ — b) . Observe that Rk belongs to M[J^i, . . . , Xn] and that Bn-p+i is well 
defined. 

Lemma 6 

For any I < ji < - ■ ■ < jp < n and I < k < n—p the polynomials (Fi, . . . , Fp, m^ '"'"'''' , 

..., M^''^ '■■■'■''' ) generate the trivial ideal or form a reduced regular sequence in 

R[Xi, . . . ,Xn]R^ ■ 

Proof. Fix 1 < ji < ■ ■ ■ < jp < n. By Condition D it suffices to prove that for 
any 1 < A; < n - p the polynomials (Fi, . . . , Fp, m[^''-'^''\ ..., mI^''-'^"^) define in 
C2^. . ./J an equidimensional variety of dimension n—p— A; and, that each irreducible 

component of this variety contains an (Fi, . . . , Fp, ]\/[{'^'"'''''' , . . . , M^''^''"'"''') -regular 
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point. The first statement follows from Lemma H] and O whereas the second one is a 
direct consequence of Conditions B and C □ 

To indices 1 < ji < ■ ■ ■ < jn "^ n we may associate a Hessian matrix Hj^j^ of 
G on Va whose entries belong to ffi[Xi, . . . ,X„]a • Let x be an arbitrary 

point of Vaj^ j • One verifies easily that the conditions det i/ji,...jp(x) 7^ and 

det J(Fi, . . . , Fp, m[^^-'^''\ ..., M£p-'^''^)(x) ^ are equivalent. 

Lemma 7 

Let I < ji < ■ ■ ■ < jp < n. The set of isolated local minimal points of G on (Vm)ajj,.,.,j 
is exactly the set of points of (Wr)a^j,..,j where i^ji,...,jp is positive definite. 

Proof. From Lemma [5] one deduces easily that the points of (VFr)a where 

^ji,--;jp is positive definite, are local minimizers of G on (Vk)a • So, we have 

only to show that the isolated local minimal points of G in (li)A belong to 

(W^k)a and that their Hessians are positive definite. 

Let X G (MR)Ai , be an isolated minimal point of G in (1^)a, , • Then, as we 
have seen, x belongs to (Wr)ajj,...j- • Let G be an arbitrary irreducible component of 
Wa ■ which contains x . Let n — p — k with l<k<n — phe the dimension of C . 

Suppose that 1 < k < n — p holds. Conditions B and C imply now that there exists 
an open subset of (Fi, . . . , Fp, M{'^'"''^''' , . . . , M^-'^''"'-^^^) -regular points of G^ which 
is dense in C^ . This implies that any neighborhood of x in Ck contains a point y 
different from x . Since G{y) = G{x) holds by Lemma the local minimal point x 
of G in (Mr) A cannot be isolated. Therefore, we have k = n — p. From Con- 

dition B we deduce that X is (Fi,...,Fp,M{^'''-'^'''\ . . . , M^^lp-'^" ^regular. Hence, 
det ifj^,...jp(x) 7^ 0. The Morse Lemma (see [H]) implies now that Hji,...,jp{^) must be 
positive definite for x being an isolated local minimal point of G in (Mr) a • ^ 

3.1.1 Unconstrained optimization 

We illustrate our argumentation in the case of unconstrained optimization. In this case 
we have p := and V is the complex afiine space C" . There is given a polynomial 
G e M[-^i, . . . ,Xn\ and the task is to characterize the local minimal points of G in 
M" . Such a local minimal point belongs to PM := {^- = 0, . . . , ^- = 0} . For the 
unconstrained optimization problem we consider the following two conditions for any 
1 < A; < n. 

Let Dk be the union of all irreducible components of W of dimension strictly larger 
than n — k . 



Condition E 



Any point 0/ { J^ = 0, . . . , j|: = 0}m \ (Dk)^ is (J^, • • • , J^) -regular. 

Let Ci, . . . ,Cs be the irreducible components of {^- = 0, . . . , ^- = 0} \ {D^) 

Condition F 

For any 1 < j < s the semialgebraic set (Cj)r \ (Ci U ■ ■ ■ U Gj-i U Cj+i U ■ ■ ■ U Gs)i 
is dense in {Gj)^. and {Gj)m. is non empty. 
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If these conditions are satisfied by G we can prove in the same way as in case of Lemma 
m [5l |6] and [7] the following corresponding statements. 

Lemma 8 

Let notations be as before. For any 1 < j < s , we have dim Cj = dim(Cj)]R = n — k . 

Lemma 9 

Let C be an irreducible component of W . Then G takes a constant value on C . 

Let 1 < k < n. By Lemma [9] the polynomial G takes on Dk only finitely many values 
which are all real. Denote by B^ the set of these values and let Rk '■= YlbeB (^ ~ ^) ^ 

M[Xi, . . . ,Xn] . 

Lemma 10 

For any I < k < n, the polynomials ^-, . . . , ^- form a reduced regular sequence in 

]R[Xi, . . . ,X„]rj. . 

Lemma 11 

The set of isolated local minimal points of G in M" is exactly the set of points of W^ 
where the Hessian of G is positive definite. 

3.2 Algorithms 

Let notations and assumptions be as in the previous subsection. We associate with 
G and Fi,...,Fp intrinsic invariants that control the complexity of the algorithms 
we are going to develop in order to solve the computational problems of minimizing 
locally and globally G on the set of points in M" defined by the equality constraints 
Fi=0,...,Fp = 0. 

Let G and Fi, . . . , Fp E Q[^i, . . . , X„] be given as outputs of an essentially division- 
free arithmetic circuit /3 in Q having size L. Let t/ > 2 be an upper bound for 
degG',degFi,...,degFp. 

3.2.1 Minimal point searching problems 

We first consider the problem of finding all isolated local minimal points of G in 
Vk = {Fi = 0, . . . ,Fp = 0}k. For this purpose we search for every index sequence 
1 < ji < ■ ■ ■ < jp < n the isolated local minimal points of G in the corresponding 
chart (1^)a • Let 5j^,,,j^ be the maximal degree of the Zariski closures in C" of 

all locally closed sets 



{Fi = 0, . . . , Fj = 0}A,^,..„,^.det//,^,....,p , 1<J<P, and 
{F^ = 0,...,Fp = 0, m[^'-'^'^ = 0, . . . , M^^'^'-*) = 0}A,,....,^.detH,,...,,, 1 < A; < n-p, 



where Hj, ~ is the Hessian matrix of G on Va 
Let finally 

5 := max{Sj,^...j^ \ I < ji < ■ ■ ■ < jp < n}. 
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We call d the degree of the isolated minimum searching problem for G on V^ = {Fi = 
0, . . . , Fp = 0}iR . From the Bezout Inequality we deduce 

Fix now I < ji < ■ ■ ■ < jp < n and observe that VFaj^,,..,j ■detHj^^,,,j is empty or a 
zero-dimensional variety. Hence, for l<k<n— pno irreducible component of the 
variety 

{Fi = 0, . . . , F, = 0, Mi^-*) = 0, . . . , m(^-*) = 0}A,,,...,,^.detH,,,...„, 
is contained in W/^. -dctH ■ In the same way as in the proof of Lemma Owe 

conclude now that Fi, . . . , Fp, MJ"'^''"'"''' , . . . , M^''^''"'''' generate the trivial ideal or 
form a reduced regular sequence in ]R[Xi, . . . , X„]a dctH • 

Applying now the Kronecker algorithm [20] to the input system 

Fi = 0, . . . , Fp = 0, Mi'"-''"^ = 0, . . . , Mi'l;-''"^ = 0, A,„...,,^ ■ det %„,.,,^ ^ 

represented by the circuit /3 we obtain for the at most 5ji,...jp complex points of 

{Fi = 0,...,Fp = 0, Mi'"-''^^ = 0,...,M£p-*^ = 0}A,,,...,,^.detH,,,...„, an algebraic 
description over Q. For the real points among them we obtain even a description a 
la Thom. We now discard the points with non-zero imaginary part and evaluate the 
signature of the Hessian matrix Hj^^j^ at each of the real points and discard the real 
points where the Hessian is not positive definite. The remaining real points are by 
Lemma [7| exactly the isolated local minimal points of G in (V]r)a • Repeating 

this procedure for each 1 < ji < ■ ■ ■ < jp < n we obtain all isolated local minimal 
points of G in ViR • 

The complexity analysis of [20] yields a time bound of L (") [n d)"^^^^ 5^ = (n d)'~'^"'^ for 
the whole procedure. Combining the Binet-Cauchy formula [17] with a similar argu- 
mentation as in the proof of [2H] Theorem 4.4, one can show that one may find proba- 
bilistically regular matrices Ai, . . . , An-p+i G Z'"^" of logarithmic heights 0{nlogdn) 
and p-minors Ai, . . . , A.„_p+i of J(Fi, . . . , Fp) ■ Ai, . . . , J(Fi, . . . , Fp) ■ An-p+i such 
that Vai U ■ ■ ■ U Va„_p+i is the regular locus of V (see [6] for details). Thus we may 
improve the sequential bound above to L (nd)'-^^^^ 5'^ . 

Theorem 12 

Let n, d, p, L, £, 6 G N with 1 < p < n be arbitrary and let G,Fi, . . . , Fg G 
Q[Xi, . . . , Xn] be polynomials of degree at most d satisfying Conditions B and C and 
having isolated local minimum searching degree at most 6 . Suppose that G, Fi, . . . , F^ 
are given as outputs of an essentially division-free circuit (5 in Q[Xi, . . . , X„] of size 
L. 

There exists a uniform bounded error probabilistic algorithm B over Q which com- 
putes from the input (3 in time L {nd)'^^^^6'^ = (nd)^^"'^ all isolated local minimal 
points of G in V^ = {Fi = 0, . . . , Fp = 0}k . 

For any n, d, p, L, i, 6 & N with 1 < p < n the probabilistic algorithm B may 
be realized by an algebraic computation tree of depth L {n d)'~'^^^ d"^ = [nd)'^^'^^ that 
depends on parameters which may be chosen randomly. 
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Let us now consider the problem of searching for the isolated local minimal points in the 
case of unconstrained optimization. Let H be the Hessian matrix of G . Observe that 
Conditions E and F imply that the polynomials ^- , . . . , ^- generate the trivial ideal 
or form a reduced regular sequence in Q[Xi, . . . , Xn\detH ■ Let 6 be the maximal degree 
of the Zariski closure in C"" of the locally closed sets {^- = 0, . . . , ^- = 0}det_ff , 
1 < k < n. We call 6 the degree of the isolated minimum searching problem for G 
on W' . The Bezout inequality implies d < {d — 1)" < d". Applying the Kronecker 
algorithm to the input system 

9G dG , ,, . 

0, . . . , — — = 0, detH ^0 



dXi ' ■ ■ ■ ' dXr, 

we obtain the analogous statement of Theorem [T2] for the isolated minimum searching 
problem in the unconstrained case. 

The aim of the next algorithm is to compute a real algebraic point which is a global 
minimal point of G in Vj^ = {Fi = 0, . . . , Fp = 0}m if there exists one. For this 
purpose we construct for each 1 < ji < ■ ■ ■ < jp < n and each 1 < k < n — p suitable 
real algebraic sample points of 

(Vi)A,, ,, n {Mi'''-''^^ = 0, . . . , M^^'^'-*) = 0}m \ {Dk)u. 

Let a; be a minimal point of G in V^ and let a := G{x) . Then x belongs to W and 
therefore there exists by Lemma [5] an irreducible component of W where G takes only 
the value a . This fact will guarantee that we are able to find a minimum oi G on V^- 
For 1 < ji < ■ ■ ■ < jp < n and 1 < k < n — p let 5ji,...jp;fc be the maximal degree of 
the Zariski closures in C" of all locally closed sets 

{Fi = 0,...,F, = 0}a, ,„./?,, 1<J<P, and 



and all generic dual polar varieties of 



{F^ = 0,...,Fp = 0, m{^'^'-*) = 0, . . . , MJ^,''-''^^ = 0}a„,...„.r„ I < k' < k 



{Fi = 0, . . . , F, = 0, Mi''-'^^ = 0, . . . , Mi''-'^^ = 0}«, 

(take Lemma [6] into account). Let finally 

S := max{5j,^,„j^.k l'^ < ji < ■ ■ ■ < jp < n, I < k < n - p} 

We call 6 the degree of the global minimum searching problem for G on V^. From 
the Bezout inequality we deduce 

6< (n(p + l)rf)" = {ndf^''\ 

We construct now recursively in l<k<n — pan ascending chain of finite sets Y^ of 
real algebraic points of Vu such that G{Yk) contains the set B^+i (see Subsection 13. II 
for the definition of -B^+i )• 
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In order to construct Yi we apply the algorithm [21 13] to the system Fi = ■ ■ ■ = Fp = . 
The algorithm returns a finite set Yi of algebraic points of V^ . 

Let now 2 < k < n — p and suppose that we have already constructed lfc_i subject 

to the condition B^ C G(Y'fc-i) • 

Let Rk := Yly^^y i^^^iv)) ^^'^ observe that Rk is a multiple of Rk in M[Xi, . . . , Xn] ■ 

Taking into account Lemma [6] and that Rk is a multiple of Rk we apply now the same 
algorithm to the system 

Fi = 0, . . . , Fp = 0, mI""-''"^ = 0, . . . , mI'''-''^^ = 0,Rk^ 0. 

In this way we obtain finitely many algebraic points of V^ which together with lfc_i 
form Yk . 

Let us consider an arbitrary irreducible component C oi W oi dimension n — p — k 
on which, by virtue of Lemma [5l the constant value of G does not belong to G(Yfc-i) • 
Then we have C^ = C and therefore C and Cr are contained in 

{Fi = 0, . . . , Fp = 0, mJ^'^'-*^ = 0, . . . , M^^'^'-*^ = 0}j^^ 

Hence, by [4], Theorem 1 the generic dual polar varieties of this locally closed set are 
not empty and the {n — p — k)th one contains a point of Cm . Applying this argument 
inductively we see that G{W) C G{Yn-p) holds. We suppose now that G reaches a 
global minimum on {Fi = 0, . . . , Fp = 0}ir . Then l^_p must contain a global minimal 
point of G in 1^ = {Fi = 0, . . . , Fp = 0}ir which is an element, say y , of l^_p with 
G{y) = min^gy„_^ G{x) . 

The complexity analysis of the algorithm of [21 [3] yields a time bound of L('^)2"~^(n d)'^^^^'^ 
(nd)'^^"^ for the whole procedure. We may use the same argumentation as before The- 
orem [T2] in order to improve this bound to L (nd)'^^^' S^ . We obtain now the following 
complexity result. 

Theorem 13 

Let n, d, p, L, i, 6 G N with 1 < p < n be arbitrary and let G, Fi, . . . ,Fs G 
Q[Xi, . . . ,Xn] be polynomials of degree at most d satisfying Conditions B, C and 
D and having global minimum searching degree at most 6 . Suppose that G,Fi, . . . , Fg 
are given as outputs of an essentially division- free circuit /3 in Q[Xi, . . . ,Xn] of size 
L. 

There exists a uniform bounded error probabilistic algorithm C over Q which computes 
from the input f3 in time L (nd)^^^^ S'^ = (nd)^^"^ a global minimal point of G in 
Vr = {Fi = 0, . . . , Fp = 0}k if there exists one. 

For any n, d, p, L, £, 6 G N with 1 < p < n the probabilistic algorithm C may 
be realized by an algebraic computation tree of depth L (nd)'^^^^6^ = (nd)'^^^^ that 
depends on parameters which may be chosen randomly. 

Mutatis mutandis the same statement is true for the unconstrained optimization prob- 
lem. The degree of the minimum searching problem in this case is the maximal degree 
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of the Zariski closures in C" of all locally closed sets {^-, • • • , '§x~}Rki ^ < k < n 
and all generic dual polar varieties of them. 

The reader should observe that Theorem [13] does not answer the question whether G 
reaches a global minimum on {Fi = 0, . . . , Fp = OJk and can only be apphed when 
this existence problem is already solved. For this problem we refer to [22]. 

4 Conclusion 

Together with [31] this paper represents only a first attempt to introduce the viewpoint 
of intrinsic quasi-polynomial complexity to the field of polynomial optimization. For 
this purpose we used some restrictive conditions which allow us to apply our algorithmic 
tools. In the future we shall relax these restrictions and extend the algorithmic tools 
and the list of real problems which can be treated in this way. 
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